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CN ! ABSTRACT 

We present a model for the satellites of the Milky Way in which galaxy formation 
is followed using semi-analytic techniques applied to the six high-resolution N-body 
simulations of galactic halos of the Aquarius project. The model, calculated using 
the GALFORM code, incorporates improved treatments of the relevant physics in the 
ACDM cosmogony, particularly a self-consistent calculation of reionization by UV 
C*~) . photons emitted by the forming galaxy population, including the progenitors of the 

central galaxy. Along the merger tree of each halo, the model calculates gas cooling 
(by Compton scattering off cosmic microwave background photons, molecular hydro- 
gen and atomic processes), gas heating (from hydrogen photoionization and supernova 
energy), star formation and evolution. The evolution of the intergalactic medium is 
followed simultaneously with that of the galaxies. Star formation in the more massive 
progenitor subhalos is suppressed primarily by supernova feedback, while for smaller 
^ \ subhalos it is suppressed primarily by photoionization due to external and internal 

sources. The model is constrained to match a wide range of properties of the present 
day galaxy population as a whole, but at high redshift it requires an escape fraction of 
UV photons near unity in order completely to reionize the universe by redshift z>8. In 
the most successful model the local sources photoionize the pre-galactic region com- 
pletely by z ~ 10. In addition to the luminosity function of Milky Way satellites, the 
model matches their observed luminosity-metallicity relation, their radial distribution 
and the inferred values of the mass within 300 pc, which in the models increase slowly 
but significantly with luminosity. There is a large variation in satellite properties from 
halo to halo, with the luminosity function, for example, varying by a factor of ~ 2 
among the six simulations. 
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1 INTRODUCTION 



E-mail: afontOast.cam.ac.uk 



A basic prediction of the cold dark matter (CDM) theory of 
structure formation is that galactic dark matter halos grow 
by the accretion and disruption of smaller subsystems. A 



2 A. S. Font et al. 




well-known consequence of this is that the Milky Way to- 
day should be surrounded by thousands of small subhalos, 
in apparent contradiction with the much smaller number 
of luminous satellites that have been detected around the 
Milky Wa y so far. This has been termed the 'miss ing satellite 
problem' ijKlypin et al] 1 19991 : iMoore et al.lll999h . A second 
potenti al discrepancy with t he theory is the recent observa- 
tion bv lStrigari et all l|200Sh that the satellites of the Milky 
Way all have very similar central densities even though they 
span five orders of magnitude in luminosity, suggesting the 
existence of a preferred scale which is not present, for exam- 
ple, in the primordial cold dark matter power spectrum of 
density perturbations. 

It has long been recognized that the reionization of 
the intergalactic medium (IGM) by the metagalactic UV 
radiation field at early times can i nhibit the formation of 
small galaxies (ICouchman fc Reed Il98fj : lEfstathioul 1 19921 : 
iThoul fc Weinberg 1996T ). It has also been recognized for 
some time that this process would provide at least part of the 
solut io n to the missing satellite probl em (Kauffman n et all 
19931: iBuTlock. Kraytsov fc Weinberg! l2000l: iBenson et al. 

Somervilld 120021 ). iBullock. Kravtsov fc Weinberg] 
in particular, emphasized the vital role of early reion- 
ization and calculated a simple model for the abundance 
of satellites in the CDM cosmogony. Since reionization im- 
poses a scale in the problem — a minimum entropy for the 
reionized hydrogen — this process may also solve the second 
'satellite problem' luTet al.ll2009l : lOkamoto fc Frenk lliooj 
IStringer et ai1l201fj| ). 

Galaxy formation in small halos can also be strongly 
inhibited by winds generated by the inje ction of supernova 
energy into the gas ( White fc Rees Ifl978l ) . The combined ef- 
fects of this so-called " supernova feedbac k " and early reion- 
ization were used by iKauffmann et all l)l993h to explain 
the relative paucity of satellite galax ies around th e Milky 
Way, and modelled in greater detail bv lBenson et al.l |2002b) 
who showed that this relatively simple physics could ac- 
count not only for the abundance of the 'classical satel- 
lites', but also for many of the properties known at the 
time, suc h as gas content, met allicity and star formation 
rate. The Ben son et al.l (|2002bh model predicted the exis- 
tence of many more satellites in the Local Group. A few 
years later, a new population of satellites was discovered, 
through careful searches of the Sloan Digital Sky Survey, 
more than doubling the number of known satellites around 
the Milky Way (IWillman et al.l 120051: IZucker et all | 200E : 
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Bclokurov et al. 20071: Irwin et al.l 120071: IWalsh et ail 12007 



see also lMartin et al.ll2004r ). The luminosity function of the 
newly found sat ellites matches the t heoretical predictions 
remarkably well (|Koposov et al.ll2008r ). 

The discovery of a new population of ultra-faint satel- 
lites led to renewed interest in the physics of dwarf 
galaxy formation. A number of recent studies have re- 
visited the arguments of the 1990s and ea rly 2000s, ei- 
ther using semi-analytic modelling techniques dCooper et all 
l2010l : iGuo et al.ll2011al: iLi et al.ll2oTol: IMaccio et al.ll2010l) 



simpler models ( Koposov et al.ll2009l : iBusha et al I l2010l : 
iMufloz et alj I2009T I . generally confirming the conclusions 
of the earlier work. Some of these studies |Busha et al.l 
l2010l : iMufioz et all 120091 ; ICooper et aT] |2010| ) have taken 
advantage of a new generation of high resolution N- 
body simulations which track the formation of galactic 



halos and their surviving subhalos down to the small- 
est halos likely to support the formation of faint galaxie s 
l|Diemand. Kuhlen fc Madaul 120071 : ISpringel et ail l2008bh . 
This allows calculation of the expected radial distribution 
of satellites which was not predicte d with sufficient prec i- 
sion by the semi-analytic approach of [Benson et al.l(|2002bh . 
based on Monte-Carlo halo merger trees and an analytic 
model of satellite orbital evolution. In addition to these 
semi-analytic calculations, full N-body/gasdynamic simu- 
lations of gala xies and their sat e llites have now been car- 
ried out fcibeskind et all 120071: lOkamoto fc Frenk 1 120091 : 
lOkamoto et al. I l20ld : IWadepuhl fc Springell l201ll v These 
too find that reionization plays a key ro le in suppressing 
dwarf galaxy formation but, as stressed by Kauffmann et al.l 
ill993h. feedback from supernova energy is also a major fac- 



Wadepuhl fc Springell l|201ll ) find, further, that cosmic 
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ray pressure is also important in suppressing star formation 
in these galaxies. 

Although a consensus seems to be emerging then, 
that the abundance and other properties of the Milky 
Way satellites can be understood as a consequence of the 
known physics of galaxy formation in a ACDM universe, 
a number of important unce rtainties remain. For example, 
iBovlan-Kolchin et alj (|201lh have recently concluded, on 
the basis of dynamical data, that the most massive sub- 
halos in N-body simulations are too concentrated to be able 
to host the brightest satellites of the Milky Way. There 
are also significant modelling uncertain ties. Recent stud- 



ies treat reionization in a simplified wa y (Busha et al 
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Guo et al.ll2011al : lLi et al.ll2010l ). reducing this complex pro- 



cess to a simple rule involving a few adjustable parameters 
(e.g., a single redshift for reionization and a threshold mass 
above which reionization is unimportant). Furthermore, the 
luminosity function is just one of many observables that de- 
pend critically on the assumed physics of the problem. The 
metallicity of the satellites, for example, not considered in 
most of these studies, is very sensitive to the effects of su- 
pernova feedback. 

In this paper, we investigate the formation of satel- 
lite galaxies in the ACDM cosmogony emp loying a new 
treatm ent of reionization, based on that of IBenson et al.l 
(|2002al ). that simultaneously solves for the coupled evolu- 
tion of the galaxy population and the IGM. The ionizing 
background is computed from the integrated emission his- 
tory of all galaxies in the Universe plus a simple model for 
the quasar contribution, filtered through the optical depth 
of the IGM whose evolution and ionization state are com- 
puted self-consistently. The ionizing radiation, generated by 
galaxies and quasars, feeds back upon subsequent galaxy 
formation through the dual action of photoionization, which 
1) prevents the collapse of gas onto low mass dark matter ha- 
los and ii) reduces the rate of radiative cooling inside more 
massive halos. An original feature of this model is the inclu- 
sion of internal reionization from progenitors of the Milky 
Way. 

Since the early work on satellite galaxies a decade ago, 
there have been important developments in techniques for 
modelling galaxy formation. Many of these have been incor- 
porated into the Durham semi-anal ytic model, GALFORM 
l|Cole et al.l 120001 : iBower et alll2006r ). which we use in this 
paper. To this, we add treatments of cooling and heating 
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processes that are relevant on the scale of dwarf galaxies. 
In addition to the usual atomic radiative processes, we in- 
clude cooling due to molecular hydrogen (H2) and Compton 
cooling of gas due to scattering off cosmic microwave back- 
ground photons. Hea ting due to photoio nization is calcu- 
lated using CLOUDY (|Ferland et alj|l998t ). The evolution of 
the IGM is computed taking into account the contribution of 
the UV ionizing background from both galaxies and quasars. 
The emissivity from galaxies is calculated self-consistently, 
whereas the contribution from quasars is modelled using an 
updated version o f the s pectrum inferred observationally by 
lHaardt fc Madaul (|l996T ). The model then predicts the time 
dependence of the properties of the IGM, such as its mean 
temperature and ionization state which, in turn, determine 
the rate at which it cools into galaxies. 

We find that, in addition to the global UV background 
pervading the Universe, the radiation generated locally by 
stars forming in the progenitors of the Milky Way is an im- 
portant source of ionization in the Milky Way at early times. 
As we will demonstrate, this radiation increases the redshift 
of reionization locally and leads to a stronger suppression 
of sate llite galaxy formatio n than the metagalactic flux. Re- 
cently, iMufioz et alj (|2009T ) have identified this mechanism 
of 'inside out reionization' of the Milky Way as important 
in suppressing the formation of ultra-faint dwarfs. Accord- 
ing to these authors, self-reionization occurs much before 
the time when the Universe as a whole was reionized, at an 
epoch when most of the gas was in molecular (-H2) form. 
They suggest that reionization would quench cooling of the 
H2 gas, stopping it from fragmenting and forming stars, pos- 
sibly explaining the paucity of ultra-faint dwarfs today. In 
practice, they model this process by making similar assump- 
tions to those made in simple models of the effects of global 
reionization, namely by assuming that H2 cooling becomes 
ineffective below a given halo mass after a given redshift. 
By contrast, we implement the self-reionization of the Milky 
Way by explicitly calculating the UV flux produced by the 
progenitors of the Milky Way as a function of time. This 
takes away any freedom in the choice of redshift and halo 
mass for which this process is important, while the same is 
achieved concomitantly for the global reionization with the 
self-consistent approach. 

In this paper, we implement our model of the coupled 
evolution of galaxies and the IGM in merger trees con- 
structed from the high resolut ion simulations of g alactic 
halos in the Aquarius project l|Springel et all 1200 8a:) . The 
resolution of these simulations (with particle masses rang- 
ing from 6.4 x 10 3 to 1.4 x 10 4 M Q ) is sufficient to enable 
predictions to be made for satellites as faint as even the 
faintest dwarfs detected recently around the Milky Way. Us- 
ing the six Aquarius simulations, we can investigate halo-to- 
halo variations in the properties of the satellites. 

Our method for simulating the joint evolution of galax- 
ies and the IGM, as well as the new physics we have im- 
plemented in GALFORM, are explained in Section [2] In Sec- 
tion[3l we demonstrate how the metallicity— luminosity rela- 
tion of satellites can be used to break the model degeneracy 
between the effects of supernova feedback and reionization. 
We also explore how self-reionization changes the predicted 
abundance of satellites. We find, however, that it is not pos- 
sible for the model simultaneously to match the luminos- 
ity function and the metallicity— luminosity relation with 



the default supernova feedback prescription used in previ- 
ous implementations of GALFORM, in which the efficiency of 
the feedback increases with decreasing halo size as a power- 
law. Instead, the satellite data require that this feedback 
efficiency should saturate for halos with circular velocity, 
«circ 65 km s . In Section 4, we test this model further 
by considering predictions for the radial distribution of satel- 
lites and for the mass within 300 pc. Finally, in Section 5, we 
discuss ou r conclus i ons an d carry out a comparison with the 
results bv lLi et alj (|2010l ) using similar semi-analytic tech- 
niques. A detailed comparison with other related models is 
presented in Appendix B. 

The cosmological parameters adopted in this study are 
the same as those u sed for the Aquarius and Millennium 
l|Springel et al.ll2005l ) simulations and correspond to a flat 
ACDM cosmology with matter density, fi m = 0.25, cos- 
mological constant term, Qa = 0.75, Hubble parameter 
h — Ho/ (100 km s _1 Mpc -1 ) = 0.73, power spectrum nor- 
malization, as = 0.9, and spectral index n s = ]Q. 



2 THEORETICAL MODELS 

We begin by describing the Aquarius N-body simulations 
used to build the merger trees employed by the semi- 
analytic galaxy formation model, GALFORM. We then de- 
scribe briefly the updates that we have made to GALFORM 
itself to include physical processes that are relevant on the 
scale of low mass galaxies. 

2.1 The Aquarius halos 

In the Aquarius project, carried out by the Virgo Consor- 
tium, six galactic dark matter halos with masses comparable 
to that of the Milky Way were sim ulated at varying levels 
of resolution (ISpringel et al.ll2008al lbl). Here, we use all six 
halos at the second h ighest level of resolution ('Level 2' in 
ISpringel et al]|2008bl ff We define as 'satellites' all sub-halos 
found by SUBFIND (|Springel et al.ll200lh with more than 20 
bound particles. Merger trees were c onstructed f o r each sub- 
halo using the met hods laid out in iHellv et all <|2003fi and 
lHarker et all <|2006h . 

Table [1] provides a summary o f the relevant p r opertie s 
of the Aquarius halos, taken from ISpringel et al.l (|2008bh . 
The halos all have virial masses in t he range 1-2 x 10 1 2 M Q . 
Recent evidence from the analysis of lGuo et alj (l2011al ) sug- 
gests that the current best estimate of the Milky Way's halo 
mass is 2 x 10 12 M Q , with a 10% to 90% confidence inter- 
val spanning the range 0.8 to 4.7 x 1O 12 M0; other studies, 

1 These parameters were chosen to match the WMAP lst-year 
results. Of these, only its and n s are of importance for our study. 
The adopted value of erg leads to an earlier redshift of reionization 
than expected for the WMAP 7-year value (erg = 0.80), but this 
is partly compensated for by the larger value assumed for n s 
(0.96 from the WMAP 7-y ear data) . In any case, these differences 
are ex pected to be small (7Boylan-Kolchin_et_al ]|2010l :l Iliev et alj 
2011), especially compared with those arising from uncertainties 
in the modelling of reionization (e.g., the assumed escape fraction 
or the clumping factor). 

2 The highest level of resolution, 'Level 1', is available only for 
one of the halos. 
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Halo 


m v [Mq] 


N hr 


M 200 [Mq] 


f200 [kpc] 


M*, to t[M ] 


name 










fbk:sat/rei:G+L (/ csc = 1) 


AqA 


1.370 x 10 4 


531,570,000 


1.842 x 10 12 


245.88 


1.085 X 10 11 


AqB 


6.447 x 10 3 


658,815,010 


8.194 x 10 11 


187.70 


9.179 x 10 9 


AqC 


1.399 x 10 4 


612,602,795 


1.774 x 10 12 


242.82 


1.050 x 10 11 


AqD 


1.397 x 10 4 


391,881,102 


1.774 x 10 12 


242.85 


6.982 x 10 10 


AqE 


9.593 x 10 3 


465,905,916 


1.185 x 10 12 


212.28 


2.217 x 10 10 


AqF 


6.776 x 10 3 


414,336,000 


1.135 x 10 12 


209.21 


6.555 x 10 9 



Table 1. Properties of the Aquarius halos, as given bv lSpringel et al (2008b). Listed here are: the particle mass, m p ; the total number 
of particles in the high resolution region, JV/, r ; the virial mass of the halo, M200; and its virial radius, T2oo- The last column gives the 
total stellar mass for each Aquarius halo in the fiducial model, fbk:sat/rei:G+L (/ eS c = 1), that will be described later in the paper. The 
stellar masses in this model show a large scatter around the value estimated for the Milky Way, 4.85 — 5.5 X 1O 1O M0 by iFlvnn et al] 
j2006h . 



however, favour a mass closer to 10 12 M q (:Battagli a et al.l 
l2006al : ISmith et al]|2007l : IXue et al.ll200ct l. Thus, the Aquar- 
ius halos sample the lower end of the currently plausible 
range of Milky Way halo masses and could underestimate 
the true mass by a factor of 2-3. This would have impor- 
tant consequences for the properties of subhalos and, there- 
fore, for any satellite galaxies they may contain. We have 
checked whether the central galaxies in our models have stel- 
lar masses comparab le to those estimat ed for the Milky Way, 
4.85-5.5 x 10 10 M Q (|Flvnn et al.ll2006T l. The final column of 
Table Q] shows the new fiducial model that will be described 
later in more detail. The total stellar masses seem realistic, 
although there is significant variation from halo to halo. 

2.2 GALFORM: The basic model 

Our starting point is the implementation of GALFORM by 
iBower et"aiT (120061'). GA L FORM includes the treatments, in- 
troduced by ICole et all (|2000h . of shock-heating, radiative 
cooling of gas within dark matter halos (modulo the reduc- 
tion in cooling rates due to photoionization described be- 
low), star formation, spheroid formation (through both disk 
instabilities and galaxy merging), chemical evolution and 
dust extinction. In addition, this vers ion contains prescrip- 
tions introduced bv l Bower et al.l |2006t ) for feedback from ac- 
tive galactic nuclei (AGN) and supernovae explosio ns, and 
mass loss in stellar winds. The I Bower et al.l (|2006t ) galaxy 
formation model reproduces properties of the galaxy popu- 
lation over a wide range of scales and epochs. For example, it 
gives a good match to the global luminosity and stellar mass 
functions and to the bimodal colour distribution of galaxies 
observed in the SDSS and can also account for the redshift 
evol ution of these pro perties . 

iFont et al] (|2008l 'l have found that the relative distri- 
bution of SDSS galaxy colours amongst the blue and red 
sequences, as well as the zero-point of the bimodal colour 
distribution, can be better matched by GALFORM by set- 
ting the m etallicity yiel d to p — 0.04 (i.e., twice the value 
adopted bv lBower etail (2006)). In addition, increasing the 
yield by this factor results in a bet ter match to the m etal- 
licity of the intracluster medium (|Bower et al.l 120081 ') and 
improves the predicted metallicities of the dwarf galaxies 
that are of interest here. For these reasons, we adopt this 
hi gher value o f the yield in this paper. As in the study 
of IBower et al.l (|2006h . we assume that the hot gaseous ha- 
los of satellites are completely and instantaneously stripped 



by ram pressure as soon as they c ross t he virial radius of 
their host halo. (|McCarthv et al.l (|2008l 'l have shown that 
ram pressure stripping of hot gas around satellites in groups 
and clusters can occur on a relatively long timescale but for 
satellites of galaxies like the Milky Way, the assumption of 
instantaneous stripping is a good approximation.) 

Feedback from supernovae plays an important role in 
establishing the properties of Milky Way satellites. In this 
paper, we investigate different feedback schemes and how 
they affect both the luminosity function and the luminosity 
- metallicity relation. The first feedback scheme we analyse is 
the im plementation in GALFORM presented in I Bower et al.l 
J2006I) (hereafter t he "default model" or fbk:B06). In the 
IBower et all |2006l ) model, supernovae are assumed to inject 
energy into the cold gas in the disk of the galaxy, heating it 
to the virial temperature of the halo after which the gas is 
ejected. This 'reheated' gas may subsequently cool and re- 
settle into the disk. The efficiency of stellar feedback, /3, is 
assumed to depend on the circular velocity of the halo, vh , 
in which the gal axy resides as /3 — (vh /vhot)~ atlat (see also 
ICole et al.ll2000T ). The paramete rs «hot and Phot con trol the 
feedback efficiency, and in the I Bower et al.l ( 20061 ) model, 
they are set to Vhot = 485 km s _1 and Ohot = 3.2. Note that 
larger values of «hot correspond to a greater effectiveness of 
supernova feedback in halos with vh < Vhot, which is the 
regime of interest here. 

As we shall see in Secti on 13.11 the stellar feedback as- 
sumed in IBower et all (2006) leads to model satellites with 
metallicities lower than observed. This happens because the 
supernova feedback in this model is too efficient at expelling 
metals from dwarf galaxies. We find that simply decreasing 
the slope of the &{vh) dependence is not a viable solution 
to this problem: while reasonable matches to the satellite 
metallicities can be found with, for example, cVhot — 2.5, 
this approach leads to a boost in the number of field dwarf 
galaxies which is inconsistent with data from large-scale sur- 
veys. 

We consider the possibility that supernova feedback 
may not scale as a simple power law across the full range 
of galaxy masses and investigate alternative formulations 
that preserve the good a greement with large-scale data that 
was the trademark of the lBower et al.l (|2006h model. In Sec- 
tion [3]2] we find a viable model which has the property that 
the efficiency of stellar feedback, /3, saturates in small mass 



halos, with v c 



^ 65 km s 



while above this value it 



retains the behaviour of the default model. This model is 
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ter subhalos are of 'type 1'. At the resolution of Aquarius 
(~ 1O 4 M ), we can safely assume that whenever a subhalo 
is disrupted so is the galaxy (if any) that it contains, that 
is, there are no type-2 galaxies. Changes in the structure of 
subhalos caused by baryons cooling inside them are likely 
to be negligible since satellites have very large mass-to-light 
ratios. We have checked that the sizes of potential type-2 
satellites given by GALFORM (including the effects of adi- 
abatic contraction caused by the bayrons) immediately be- 
fore their subhalos are disrupted are generally larger than 
the subhalo tidal radius. We therefore assume that 'type 2' 
galaxies are disrupted when their halos are disrupted and 
exclude them from further analysis. 

With the exception of the additional physical pro- 
cesses discussed below, all the oth er parameters in our semi- 
analytic model are the same as in iBower et al l (|2006l ). 

2.3 Additional physics 

Here, we briefly describe the inclusion of new physical pro- 
cesses, relevant on dwarf galaxies scales, in the GALFORM 
model. 



Figure 1. Comparison of the global K-band luminosity function 
for the fiducial mo del with satu r ated feedback used in this paper 
(full line) and the IBower et ay 1120061 ) model (dashed line). The 
model s are compared with observational data from ICole et al.l 
(2001). Although the two models differ in their treatment of stel- 
lar feedback below i> c irc = 65 km s _1 , they provide entirely com- 
parable fits to the observational data at brighter magnitudes. 



called fbk:sat. It matches the global properties of the galaxy 
population well (see Fig. [1] for example) and, as we will 
show below, also a variety of Local Group data, including 
the flattening of the metallicity-luminosity relation observed 
for the ultra-faint dwarfs. The stellar feedback is required to 
be very inefficient in these systems. 

It seems plausible that the feedback efficiency is low 
in small mass galaxies because their star formation rate 
(and hence the supernova rate) is observed to be low. Dwarf 
galaxies may be pre-en riched by Popula tion III stars to a 
level of order 10" 3 Z Q (jWise et all |2010| ). Afterwards, due 
to the inefficient feedback, the smallest dwarfs would evolve 
appro ximately as a closed-box system fSalvador i fc Ferraral 
120091 ). and a single generation of supernovae may be suffi- 
cient to increase their metallicity by ~ 0.5 dex, as observed. 
We note that a s imilar satu r ation scheme has been pro- 
posed recently by iGuo et al.l (|2011ah . who implemented it 
in the Munich semi-analytic model. Theirs is the first model 
that successfully matches both the large-scale data and the 
properties of dwarf galaxies, although only at the resolu- 
tion of the Millennium-II s imulation, m p ~ 9.4 x 1O 6 M0 
l|Bovlan-Kolchin et al.ll2009l ). Tabled summarizes the main 
feedback and reionization features of the models used in the 
paper. 

When modelling galaxy formation in lower resolution 
simulations than Aquarius care must be taken to follow 
satellite galaxies in subhalos that are disrupted below the 
adopted particle li mit (20 in our case) by tidal forces. In 
the terminology of lSpringel et all l|200ll ). these galaxies are 
of 'type 2' while galaxies present in identified dark mat- 



2.3.1 Cooling, reionization, and photoheating 

The cooling function, A(p, T, Z, z) (which is the net cooling 
rate of gas, obtained by summing the cooling and heating 
terms), is calculated self-consistently at each redshift us- 
ing the photoheating background predicted by our model 
of the evolution of the intergalactic medium (IGM) as de- 
scribed below. To carry out this ca lculation we employ the 
code CLOUDY (|Ferland et all 119981 ). The radiative cooling 
processes included are: Compton cooling off the cosmic mi- 
crowave background, thermal bremsstrahlung, and the usual 
atomic radiative processes. We also include molecular (H 2) 
cooling following the prescription of lBenson et al.l l|200rj ). 

Our model for the evolution of t h e IGM is essen- 
tially that described in iBenson et al.l j|2llll(i : (see also 
iBenson fc Boweil l2010l ). which follows the photoionization 
and recombination of the IGM in addition to cooling and 
heating rates thereby allowing the ionization and thermal 
state of the IGM to be predicted as a function of time. The 
photoionizing flux present at any point within an Aquar- 
ius halo is made up of two contributions, a local one due 
to sources within the halo and a global one due to sources 
in the rest of the Universe. The calculation of the local 
flux is explained in Section \2. 3. 21 To obtain the global flux, 
we run the GALFORM code on a very large set of Monte- 
Carlo merger trees generated using t he em pirical modifica- 
tion of the extended iPress- Schechterl Il974l ) formalism advo- 
cated bv lParkinson. Cole fc HeTivl (|2008l ). which gives results 
consistent with N-body simulations. This method gives the 
properties of the IGM outside the Aquarius halo, as well 
as the global photoheating background. At each timestep in 
the GALFORM calculation, we determine the mean emissiv- 
ity (as a function of wavelength) per unit volume in the Uni- 
verse by summing the contributions of galaxies and quasars. 
This emissivity is used to compute the rate at which the 
background of ionizing photons is built up. The background 
can experience absorption by neutral gas in the IGM, and 
so it is strongly suppressed prior to reionization. 
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Model 


Supernova feedback 


Global 


Local 


Escape 


name 


efficiency (/3) 


reionization 


reionization 


fraction (/esc) 


fbk:B06/rei:G 


Bower et al. (2006) 


/ 


X 


1 


fbk:B06/rei:G+L 


Bower et al. (20061 


/ 


/ 


1 


fbk:sat/rei:G 


Saturated 


/ 


X 


0.1, 0.2, 0.5, 1 


fbk:sat/rei:G+L 


Saturated 


/ 


/ 


0.1, 0.2, 0.5, 1 



Tabl e 2. Features of the models. The name of each model reflects whether it includes the default supernova feedback (fbk:B06) efficiency 
of lBower et al.l (120061 1 or the saturated feedback efficiency (fbk:sat), whether it includes global (rei:G) or local (rei:L) reionization. (Details 
of the reionization model are given in !2.30 



For galaxies, the emissivity is obtained self-consistently 
from the population of galaxies formed by GALFORM up to 
the current timestep. The efficiency with which this radia- 
tion ionizes the IGM depends on the assumed escape fraction 
of ionizing photons. While both observati onal and theoret- 
ical estimates of this fraction exist (e.g. Atek et al. 20091 ; 
iTaursen et aT1l2009l ; ISiana et al.ll2007l ; I Wise fe Chenll2009h . 



its value remains highly uncertain, particularly at high red- 
shift. We therefore chose to treat the escape fraction as 
a free parameter which we fix so as to produce a reion- 
ization history comp atible with experimental constraints 
|Komatsu et al.ll2010l ). As we will show in Section [3.21 we 
find that our new model requires a high escape fraction 
(^ 80 — 100%) in order to match the present number of 
satellite galaxies. The high escape fraction reflects the fact 
that our model produces too few ionizing photons at high 
redshifts to reionize the Universe sufficiently early. 

The escape fraction in our model is clearly unreal- 
istically large. Observational est i mates suggest value s of 
<;i - 20% (e.g.lSiana et al.[l2007l ; rVanzella et al1l2010l ; but 
see lBolton fc Haehneltll2007l l. but we are forced to adopt a 
value of ~100% to achieve reionization by a sufficiently high 
redshift. While it is unclear whether the observed values 
are representative of the galaxy population as a whole (e.g. 
if unobserved, low luminosity galaxies dominate the produc- 
tion of ionizing photons, and their escape fractions are much 
larger, then the ionizing emissivity-weighted escape fraction 
may be much larger), this is nevertheless a point of tension 
between our model and current empirical expectations. 

Other models of the formation of the Local Group satel- 
lites have not attempted to model the formation of galaxies 
and the process of reionization self-consistently as we have 
here, and so have been able either to choose an epoch of 
reionization consistent with constraints, or to put in an ion- 
izing emissivity by hand that is designed to give a reasonable 
epoch of reionization. 

For the purposes of this work, which is concerned pri- 
marily with the Local Group satellite galaxies, and not with 
the process of reionization itself, we take the pragmatic view 
that we should do whatever is necessary to obtain a realistic 
epoch of reionization since it is this that primarily impacts 
the formation and suppression of dwarf galaxies. The fact 
that this requires a large escape fraction is interesting, and 
deserves further attention, but is not surprising given that 
both observations and our model are highly uncertain in the 
high redshift regime. A better model will most likely need 
to produce more ionizing photons at high redshift thereby 
allowing the escape fraction to be reduced. Of course, a self- 
consistent calculation of the escape fraction for model galax- 



ies would provide a more rigorous w ay to approach this prob- 
lem (e.g. iFernandez fc Shullll201ll ). A more detailed study 
of reioni zation in the GALFO RM model is currently under- 
way (see iRaicevic et ai1l2010l ) and will explore this issue in 
greater depth. 

For calculating the emissivity of qu asars, we use 
the code CUBA (|Haardt fc Madaul l200ll ). which in- 
cludes an updated vers ion of the observationally inferred 
lHaardt fc Madaul (|l996l ) spectrum. This approach allows us 
to track the evolution of the IGM while simultaneously cal- 
culating the heating rate of the gas inside galaxies. 

Heating due to UV photoionization (photoheating, for 
short) not only offsets cooling losses but it also prevents low 
mass dark matter halos from accreting their full comple- 
ment of baryons. The suppression of baryonic accretion into 
halos is modelled using the accr etion mas s scale recently de - 
rived bv lOkamoto et al. I (|2008l ) (see also lHoeft et al.ll2006h 
from cosmological hydrodynamical simulations. These au- 
thors find that the baryon fraction, /b = Mh/M, depends 
on halo mass and redshift as: 



f h (M,z) = (f h ) 



1 + (2 Q/3 - 1) 



M 



M c (z) 



-3/c 



(1) 



where (/b) = fib/fio is the universal baryon fraction, M c (z) 
is a characteristic mass that will be defined below, and a — 2 
provides a good fit t o the results of the si mulations. It is 
worth noting that the lOkamoto et aL 1 (120081 ) accretion mass 
scale is significantly lower than t he "filtering mass" calcu- 
lated previously bv lGnedir] (|2000l ) on the basis of linear per- 
turbation theory. With this new method, halos that accrete 
only half of the universal baryon fraction (i.e. /b = (,fb)/2) 
have, on average, circular velocities of about 25 k m s~ at 
redshi ft z = 0, compared with 50 km s" 1 found b y lGnedinl 
(2000). In other words, the lOkamoto et al. I (|200St ) mass ac- 
cretion scale allows more satellites to escape the effec ts of the 
reioni zation than the associated filtering mass in the lGnedinl 
(2000) formalism. (The latter was used in the previous ver- 
sions of GALFORM). As we will show, this has important 
consequences for the inferred role of supernova feedback on 

relatively massive satellites. 

W e follow the method suggested by lOkamoto et al. I 
(2008) for implementing their scheme of baryon accretion, 
which involves computing the equilibrium temperature of 
gas as it accretes into the halo, r acC r. The simple expecta- 
tion is that the characteristic mass scale should be set such 
that T vir (M c [z], z) = T accl (z) where T vil (M) is the virial 
temperature of a halo of mass M at redshift z. However, 
lOkamoto et al. I l|2008h show that this approach overesti- 
mates the suppression mass scale and does not accurately 
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Figure 2. Evolution of the IGM in the default feedback jBower et al.ll2006l ') model. Left: the mean temperature of the gas in the proto- 
Local Group region in the models with and without the added local UV flux, as a function of redshift (dashed and solid blue curves 
respectively). Reionization occurs at z ~ 14 in the model with local and global UV flux and at z ~ 6 in the model with only global UV 
flux. In the first model, the temperature peaks at just below 3 X 10 4 K whereas in the second model it peaks at 2 X 10 4 K. Both models 
assume an escape fraction of 100%. The temperature of the CMB is shown with a red dotted line for reference. Right: the characteristic 
accretion temperature of gas, T a c C r, as a function of redshift. The accretion temperature peaks at z ~ 3, at a value of 8 X 10 4 K in the 
model with local UV production, and only at z ~ 1.5, at a value of 4.5 X 10 4 K, in the model with only global UV production. Note the 
different scales on the y-axis of each panel. 




Figure 3. Evolution of the mean temperature of the gas in the proto-Local Group region in the model with saturated feedback. The 
left panel corresponds to the case of global reionization only, with different colours indicating the escape fraction. Global reionization 
occurs at z ~ 7.8 for / osc = 1, z ~ 6 for / csc = 0.5, z ~ 5.5 for / csc = 0.2 and z ~ 4.5 for / csc = 0.1 respectively. T proto LQ peaks at 
~ 2 X 10 4 K for / csc = 1 and decreases slightly to ~ 1.5 X 10 4 K for / CS c = 0.1. The temperature of the CMB is shown with a red dotted 
line for reference. The right panel shows a comparison of T pro toLG m the saturated feedback model with and without local reionization, 
for /esc = 1 (dashed and solid blue lines, respectively). Local emission shifts the redshift of reionization from z ~ 7.8 to z ~ 10, without 
a significant increase in the peak of the gas temperature. 
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reproduce the redshift dependence found in their simula- 
tions. Instead, they recommend a simple model in which 
each halo accretes gas at the universal rate (i.e. f2t>/f2o times 
its total mass accretion rate) if T v j r > T accr , and accretes no 
gas if T v i r < T accr . We adopt their model in this work. The 
accretion temperature is computed self-consistently using 
the cooling function described above and the known density 
of accreting material. Note that the lOkamoto et al. I (|2008l ) 
formalism provides the accretion mass scale only after the 
epoch of reionization, when there is a significant ionizing 
background, at which point the authors suggest setting the 
critical temperature to the equilibrium temperature for the 
photoionized accreting gas. At earlier times, prior to reion- 
ization, we assume that the pre-shock temperature of gas 
as it is accreted onto a halo is given by the temperature 
to which it is adiabatically heated up from the mean IGM 
temperature, or 10 4 K, whichever is lower. This is a valid 
approximation prior to reionization as the IGM gas will be 
mostly neutral, atomic and metal-free and so there are no ef- 
ficient cooling processes below 10 4 K, and little photoheating 
due to the lack of a significant ionizing background. At later 
times, the full complement of heating and cooling processes 
are used to compute the accretion temperature. 

Our treatment of reionization can be compared to 
othe r semi-analytic/semi-numerical app roaches. In partic- 
ular, iMesinger. Furlanetto fc Cenl |201ll ') recently described 
a semi-numerical algorit hm to compute the ionization state 
of the IGM. Unlike us, IMesinger. Furlanetto fc Cenl (|201lh 
compute a position-dependent ionization fraction by re- 
course to cosmological simulations which they use to esti- 
mate densities, collapsed fractions and star formation rates 
(assuming a constant efficiency to convert collapsed frac- 
tion to star formation). In contrast, our model assumes a 
uniform reionization, but additionally solves for the thermal 
state of the IGM and has a much more detailed treatment of 
the underlying galaxy formation physics to predict the ioniz- 
ing emissivity. Further advances in modelling of reionization 
must co mbine aspects of thes e and other approaches. For 
example. iRaicevic et al.l (|2010T l compute position-dependent 
reionization by combining a detailed model of galaxy forma- 
tion with ray-tracing through a cosmological simulation of 
the density field. 

2.3.2 The local UV flux from Milky Way sources 

In addition to the global photoionizing flux originating from 
large scales (from AGN and quasars), local sources can also 
be important contributors, especially at high redshift. If 
these local sources are significant, reionization in the re- 
gion will occur earlier than in an average region of the Uni- 
verse. Thus, the entire region destined to become the Local 
Group may be encompassed within an ionized bubble ( "local 
reionization"), before many such bubbles have percolated to 
reionize the Universe as a whole ("global reionization"). One 
consequence of local reionization is to suppress further the 
formation of dwarf satellites within the region. 

A fully self-consistent treatment of local reionization re- 
quires not only accounting for the emissivity of local sources, 
but also knowledge of other (at the moment, poorly con- 
strained) parameters, such as the photon escape fraction and 
the gas clumping factor, and modelling of radiative trans- 
fer. A comprehensive treatment of this kind has not been 



performed so far. Current models that do include rad iative 
transfer (e.g. IWeinmann et al.ll2007l ; llliev et all 120111 ) have 
very simplified physical prescriptions for star formation and 
feedback. Semi-analytical models such as GALFORM, which 
have more realistic star formation and feedback prescrip- 
tions and which match a wide range of observations, by con- 
trast, do not include radiative transfer. The combination of 
the two approaches is clearly desirable for a c omplete under- 
stand ing of the role of local reionization (see IRaicevic et aU 
2010 for an example of progress in this area). Below we 
present a simple method for accounting for the local pho- 
tons. Without radiative transfer, this is likely to overesti- 
mate the contribution of local photons to the suppression 
of low-mass satellites. Nevertheless, we present the results 
here as proof of concept of how the redshift of reionization 
is increased by contribution from pre-galactic sources. 

To calculate the contribution to the ionizing flux from 
local sources, we identify the galaxy's progenitors by finding 
halos in the merger tree that lie within the virial radius of 
the final Aquarius halo at z = 0. We then calculate the 
total luminosity of ionizing photons emitted by these local 
sources at all redshifts, and compute the effective photon 
density within the Lagrangian radius of the final halo: 



where Li on (z) is the total instantaneous ionizing luminosity, 
c is the speed of light and J?La g (z) is the radius of a sphere 
with volume equal to the physical volume that contains a 
mass equal to the mass of the final halo at each redshift, 
assuming that the sphere is at the mean density. To account 
for local emission, this additional emissivity is added to the 
mean emissivity computed from the global distribution of 
galaxies. We find that this local emission is the dominant 
contribution to the net background at early times, particu- 
larly prior to reionization when the global background has 
been unable to build up as a result of the high optical depth 
in the IGM. At later times, e.g. for z < 2, the local con- 
tribution becomes entirely negligible. Accounting for the lo- 
cal emission in this way is an approximation that will be 
valid when the integrated background is dominated by re- 
cently emitted light. Prior to reionization, this is certainly 
the emitted light is rapidly absorbed by neutral gas. 

Post reionization, the local contribution to the emissivity be- 
comes small compared to the mean global emissivity and so 
can be neglected. 

Fig. [2] shows the evolution of the proto-Loca l Group 
region in the default feedback (I Bower et all l2006h model, 
adopting an escape fraction of ionizing photons of 100%. The 
left-hand panel illustrates the redshift evolution of the mean 
gas temperature in the models with and without the local 
UV flux. At high redshift, the gas cools adiabatically due 
to the expansion of the Universe. As the first galaxies form, 
they begin to photoheat the IGM resulting in a rapid rise 
in temperature at z ~ 20 when the gas temperature soon 
exceeds that of the CMB by over 2 orders of magnitude. 
In the model with both global and local photoionization, 
the proto-Local Group region is reionized at z ~ 14. At 
that point, the gas reaches its maximum temperature of just 
under 3 x 10 4 K. In the model without the contribution from 
local sources, reionization occurs s ignificantly later, at z ~ 6 
(as in the previous calculation of iBenson et al.ll2002ah and 
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the maximum temperature that the gas attains is only 2 x 
10 4 K. After reionization, the temperature of the proto-Local 
Group gas declines slowly, at a rate controlled by the balance 
between cooling and continued photoheating. 

The right-hand panel of Fig. [2] shows the characteristic 
accretion temperature of the ga s, T accr , as a function of red- 
shift (see lOkamoto et al. 112008 ). The accretion temperature 
rises rapidly at the epoch of reionization as a strong photo- 
heating background builds up, reaching a value of ~ 8 x 10 4 K 
at z ~ 3 in the model with local sources. In the model 
without local sources, the maximum temperature is approx- 
imately half as high, 4.5 x 10 4 K, and this is only reached at 
z ~ 1.5. At later times, the accretion temperature begins to 
fall as the photoheating background declines. 

The left panel of Fig.0shows the evolution of the mean 
temperature of the proto-Local Group region in the model 
with saturated feedback and global reionization only, for var- 
ious assumed escape fractions of ionizing photons. As ex- 
pected, for a lower escape fraction global reionization oc- 
curs later, but the model results in a more plausible reion- 
ization history fo r large values of the escape fraction (see 
IStark et all l20ld . for a discussion of recent observational 
data). The earliest redshift of reionization in this case is 
z ~ 7.8, for / csc = 1, and the latest is z ~ 4.5, for / csc = 0.1. 
The peak of T pro toLG is similar in the two feedback models 
(compare the solid blue lines in the left panels of Figs. [2] 
and [3}. For / csc = 1, this temperature is ~ 2 x 10 4 K. At 
a fixed escape fraction, global reionization occurs earlier in 
the model with saturat ed feedback t h an in the model with 
the default feedback of iBower et al.l |2006) (z ~ 7.8 com- 
pared with z ~ 6 for / csc = 1). This is because a weaker 
feedback efficiency allows dwarf galaxies to become brighter 
and therefore emit a larger number of ionizing photons. 

The right panel of Fig. [3]shows the effect of adding local 
reionization in the model with saturated feedback. (We illus- 
trate this for /esc = 1, but the other cases display a similar 
behaviour). In this case, local emission shifts the redshift of 
reionization from z ~ 7.8 to z ~ 10, while the peak of the 
gas temperature remains roughly the same. Below, we dis- 
cuss how different feedback and reionization models affect 
the satellite luminosity function. 



3 RESULTS 

In this section, we investigate how different feedback pro- 
cesses affect the broad properties of satellite galaxies. We 
compare our model predictions not only to the observed 
satellite luminosity functiorff], but also to the observed satel- 
lite metallicity-luminosity relation. As we shall see, this re- 
lation serves to break a degeneracy in the models between 
the effects of supernova feedback and reionization. 

We will compare our model predictions for the luminos- 
ity function with data for the Milky Way, based largely on 
the SDSS DR5 data. For this, we count all satellites in the 
model within a radius of 280 kpc, the limit to which the 
tip of the red giant branch can be detected in the SDSS. 
We compared this luminosity function to the estimate by 



3 Unless stated otherwise, we use the term 'luminosity function' 
to refer to the luminosity function of dwarf satellite galaxies 
around galaxies like the Milky Way. 
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Figure 4. Luminosity functions for the models with default feed- 
back (full lines for global reionization: fbk:B06/rei:G; dashed lines 
for global plus local reionization: fbk:B06/rei:G+L). The different 
colours correspond to the six Aqu arius halos. The estim ate of the 
Milky Way luminosity function bv lKoposov et all | |200S| ) is shown 
by the stepped dotted line down to the limit of My = — 5, below 
which the volume corrections become very uncertain (this limit 
is marked in the figure by a vertical arrow). The errorbars are 
Poissonian for the satellites in the classical (My ^ — 11) regime, 
while for the fainter systems they include volume corrections ap- 
propriate for the SDSS DR5. 



IKqposov et al.l l|2008l ) of the expected number of satellites 
in the Milky Way within this radius, which is well described 
by a power-law, 



dN/dM v = 10 x i 01 ( M v+5). 



(3) 



The number of satellites in the Milky Way estimated ac- 
cording to this calculation is shown in the luminosity func- 
tion plots below by a stepped dotted line. In the "classical 
dwarfs" regime (My ^ — 11), the error bars on the obser- 
vational estimate are assumed to be Poissonian, while in 
the interval —11 ^ Mv ^ —5, additional corrections for the 
SDSS DR5 volume are included. We do not plot the observa- 
tional estimates for magnitudes fainter than Mv — — 5 where 
the volume corrections become very uncertain. Note that in 
this cumulative plot the errors are not independent. Also 
note that the above estimate requires an uncertain assump- 
tion about the radia l distribution of Milky Way satellites 
l|Tollerud et alj|2008t ). 

For properties of the satellite population for which the 
radial incompleteness of the observational data is important 
(e.g. the radial distribution of satellites in Section [4.ip . we 
model the incompleteness explicitly within the SDSS DR5 
footprint. Firstly, we apply the magnitude de pendent thresh- 
old of detectability, R m ax{My), derived by iKoposov et all 
(2008). Then we randomly select 20% of the satellites to ac- 
count for the SDSS DR5 partial coverage of the sky. The 
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Figure 5. The Z — My relation for the fbk:B06/rei:G models (left) and for the fbk:sat/rei:G+L models (right), where Z is the average 
[Fe/H]. Both models assume / csc = 1. The different colours co rrespond to satellit es in the six Aquarius halos. The black squares show 
the observat i onal data for the Milky Way dwarf sphe roidals iNo rris et al~ll201o |). This study includes data for both cl assical dwarfs 
l|Mateolll998l : iHelmi et aDl200rj) and ultra-faint dwarfs l|Kirbv et alj|2008l : iMartin. de Jong fc Rixlliopj : INorris et al~ll2010ll. b ut only for 
those syste ms with chemically un biased star samples. The [Fe/H] error bars have been re-calculated by INorris et al. using the 

pr ecepts of |Pa Costa et al.l ll977] f§ III). The star symbol shows the peak of the metallicity distribution function in For nax measured 
by iBattadia et alJ j2006bh and iKirbv et al . (2010) using a larger sample of stars but with lower resolution spectra than INorris et al. I 
The lower resolution samples cover a larger spatial extent in Fornax and therefore may give a more accurate repre sentation of the 
average [Fe/H] by taking i nto account the contribution of metal-poor populations preferentially located at the outskirts ( Battaglia et al. 
l2006bl : iLetarte et aljlioiol) . 



model "classical" satellites (My ^ — 11) are all included, 
under the assumption that they are not affected by incom- 
pleteness. 



3.1 The default feedback model (fbk:B06) 

The combined effects of supernova feedback and reionization 
in the default feedback model, applied to the six Aquar- 
ius halos, are shown in Fig. [4] Models with global reioniza- 
tion only are labelled fbk:B06/rei:G (solid lines) and models 
which, in addition, include local reionization are labelled 
fbk:B06/rei:G+L (dashed lines). 

Within the scatter, the fbk:B06/rei:G models give a rea- 
sonable match to the luminosity function. However, none 
of the halos host galaxies as bright as the the SMC and 
the LM C. A similar result was obtained by I Benson et al.l 
|2002ri ), who found that only about 1 in 20 of their Milky 
Way-type galaxies had satellites as bright as the SMC and 
the LMC. Remarkably, this result now appears consistent 
with recent measurements of the prevalence of bright satel- 
lites in external galaxies similar to the Milky Way which 
indicate that only 11% of such hosts have one and 3.5% 
two s atellites as bright as the Magellanic Clouds ijLiu et al.l 
120111 ). Note that there is a large variation in the number 
of model satellites from one halo to another: the predicted 
satellite abundance varies by a factor of about 2, reflect- 
ing the different formation histories of the halos. This rela- 
tively large scatter highlights the danger of arriving at far- 
reaching conclusions regarding cosmology based on the sin- 



gle example of the Milky Way. Indeed, iGuo et all (|2011bh 
have recently shown, using SDSS data, that, in the mean, 
isolated primaries of comparable luminosity to the Milky 
Way contain about a factor of two fewer satellites brighter 
than M v = -14 than the Milky Way itself. 

In the models, there is a degeneracy between the ef- 
fects of supernovae feedback and reionization: both suppress 
galaxy formation in small halos. The metallicity, Z, of the 
stars and gas in a galaxy that is already assembled is not af- 
fected by photoheating. It can, however, be strongly affected 
by supernova feedback which reduces the effective yield as a 
consequence of outflows. Thus, the Z — My relation has the 
potential to break this degeneracy. 

The Z - My relation in the fbk:B06/rei:G model is 
shown in the left panel of Fig. [5] for all six Aquarius ha- 
los, and compared with observations o f Milky Way dwarf s 
with reliable metallicity measurements l|Norris et al. II20T0I '). 
The default feedback model undershoots the metallicities 
of dwarf galaxies, particularly for the less luminous sys- 
tems. (Note that the m odel already includ es the large yield, 
p = 0.04, favoured by iFont et all ^OOSh ). The scatter in 
the Z — My relation for all six Aquarius systems is rela- 
tively small, suggesting that the merger history plays only 
a secondary role in determining the slope or zero-point of 
the relation. We conclude that the supernova feedback in 
this model is too efficient at expelling metals from dwarf 
galaxies. We now analyze the saturated feedback model. 
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Figure 6. Satellite luminosity functions in the fbk:sat/rei:G and fbk:sat/rei:G+L models for the six Aquarius halos. The photon escape 
fraction is taken to be / csc = 1. The solid lines represent the models with global reionization only and the dashed lines t he models with 
both local and global reionization; different lines correspond to the six Aquarius halos. The observational estimate of iKoposov et all 
(120081 ) is shown as a stepped dotted line, as in Fig. [4] 



3.2 The saturated feedback model (fbk:sat) 

The right panel of Fig. [5] shows the Z — My relation in 
the model fbk:sat/rei:G+L. With this feedback prescription, 
the match to the observed Z — My relation is greatly im- 
proved, especially for the fainter satellites. In particular, the 
turnover in satellite metallicities is well reproduced. As dis- 
cussed in Section 12.21 a floor in the average metallicity of 
satellites can be created by a very inefficient supernova feed- 
baclfl There are two main regimes in the evolution of dwarf 
galaxies with inefficient feedback. At early times, there is 
a significant gas inflow that sustains star formation and an 
increase metallicity up to ~ W~ 2 ' 5 Zq (a single generation 
of supernovae may be sufficient to increase the gas metallic- 
ity by ~ 0.5 dex). This regime is short for the smaller dwarf 
galaxies, as their gas inflow drops quickly. From then on, the 
rate of change in metals becomes directly proportional with 
the outflow which, assuming a constant feedback efficiency, 
leads to a floor in metallicity. 

At the brightest end (i.e., Fornax) our model pre- 
dicts average metal licities that are below the data point of 
iNorris et al. I l|2010r ). However, this measurement is likely to 
be an overestimate of the actual value because the sample 
used preferentially contains stars close to the centre of the 
galaxy (which have highe r resolution spectra ), and here the 
stars are more metal-rich l|Letarte et al.ll201oh . Our data are 
in better agreement with the average Z values inferred from 



4 We could not reproduce this flattening of Z at faint My with 
a typical power-law supernova feedback, regardless of the change 
in slope or zero-point. 



lower resoluti on spectra that cover a larger extent of this 
dwarf galaxy l|Battaglia et al.ll2006bl ; iKirbv et al.ll2010t ). 

Fig. HJ] shows the luminosity functions of the six Aquar- 
ius halos in the fbk:sat/rei:G and fbk:sat/rei:G+L models, 
for an assumed photon escape fraction of / eS c = 1. This high 
value is required to match the observed luminosity function. 
As mentioned earlier, this value produces a plausible reion- 
ization history, but seems higher than indicated by some 
observational data. Although the escape fraction remains 
uncertain, such high values (which arise in GALFORM be- 
cause too few ionizing photons are produced at high red- 
shifts) is probably a shortcoming of the model and requires 
further investigation. However, for our purposes here, the 
actual value of the escape fraction is not in itself impor- 
tant. What matters is that there should be enough ionizing 
radiation to suppress the formation of small galaxies. 

Fig. [5] indicates that, unlike in the fbk:B06 case, 
global reionization alone does not suppress galaxy forma- 
tion enough in the fbk:sat model to account for the satel- 
lite luminosity function. In this case, reionization by local 
sources is also required. There are important differences in 
our two feedback models. In particular, the saturated feed- 
back model produces more satellit es of all lumi n ositie s, even 
with / osc = 1, than the default iBower et all (|2006l ) feed- 
back model (compare Figs. [3] and [5]) This is largely because 
the lower feedback efficiency in the saturated model allows 
galaxies to grow brighter. 

In summary, we have found that the model 
fbk:sat/rei:G+L matches both the global luminosity 
function (of all galaxies; see Fig. [TJ, the local luminosity 
function of Milky Way satellites (see Fig. [6} and the 
Z — My relation (see Fig. [5}. This model also produces 
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Figure 7. The effect of different physical processes on the luminosity function of satellites in the Aq-A halo, in the models fbk:B06 (left) 
and fbkisat (right). Both models shown here have / csc = 1. The purple short-long dashed line shows a model without any reionization 
or supernova feedback. The dark blue dot-long dashed line shows a model with supernova feedback only, no reionization. The two pure 
dashed lines show models without supernova feedback but with reionization, dark red long dashes for global reionization only and orange 
short dashes for global plus local reionization. The cyan full line shows the model with global reionization and feedback. The model with 
global plus local reionization and supernova feedback is shown with full dark green lines. The observational data and their errorbars are 
as in Fig. [4] 



central galaxy stellar masses roughly consistent with the 
measured stellar mass of the Milky Way (see Table HJ. In 
Section [4] we will adopt this as the fiducial model and show 
that it also matches other properties of the Milky Way 
dwarf satellites. 

3.3 The effects of supernova feedback and 
reionization 

Before exploring our fiducial model further, it is instructive 
to compare the separate effects of the two sources of feed- 
back, supernova energy and reionization, on the satellite 
luminosity function in the two cases, fbk:B06 and fbkisat. 
These are illustrated in the two panels of Fig. [7] 

In the absence of feedback of any kind, the purple dot- 
dashed lines show that in both cases hundreds of satellite 
galaxies form with ^/-magnitudes brighter than My = —10. 
The lack of feedback enables galaxies to retain their baryons 
and continue to grow in size and luminosity. Below My « 
— 10, the cumulative luminosity function levels off due to the 
inability of gas in halos with virial temperatures below 10 4 K 
to cool efficiently. The effect of reionization by a global UV 
background (dark red long dashed lines) is not enough, on 
its own, to lower the number of satellites sufficiently so as to 
match the observed luminosity function. The effect is mini- 
mal in the fbk:B06 model, where it suppresses the formation 
of satellites only by about 50%. In the fbk:sat model it re- 
duces the number of luminous satellites by a factor of 2, but 
still overpredicts the bright end. Our model of reionization is 
based on the hydrodynamic simulations of lOkamoto et al~l 
(2008) in which the effects are relatively mild, with gas ac- 



cretion being suppressed only in halos with circular velocity 
lower than ~ 25 km s - . Most previous studies of satel- 
lite galax ies have assum ed the more aggressive reionization 
model of lGnedinl £2000) in which suppression occurs in ha- 
los with circular velocity up to ~ 50 km s _1 . However, even 
in this case, the effects of global r eionization are relatively 
mild (see e.g. lBenson et alj|2002bh . 

The inclusion of local sources of reionization (orange 
short dashed lines in Fig. [7]) has a dramatic effect in fbk:B06, 
but a relatively weaker effect in the fbk:sat model. This is 
to be expected since in the fbk:sat model a larger fraction of 
ionizing emission arises from lower mass progenitors. In the 
overdense region corresponding to the proto-Local Group 
these low mass progenitors are less overabundant (relative 
to their mean abundance averaged over all s pace) than are 
more massive progenitors l|Mo fc Whites 1996). Thus, the en- 
hancement in the total ionizing emissivity is similarly re- 
duced when accounting for their local contribution. In both 
cases reionization is pushed back locally to an earlier red- 
shift and galaxy formation is inhibited in halos with very 
low velocity dispersion. The bright end, My < —15, remains 
unaffected. 

Perhaps the most illuminating is the effect changing 
the supernova feedback. In the fbk:B06 model this process 
affects galaxies in halos of all circular velocities. The ef- 
ficiency is high and supernova feedback dominates global 
reionization across the entire range of satellite galaxies. In 
the fbk:sat model, this process becomes scale-dependent. In- 
terestingly, the scale on which supernova feedback saturates, 
My w —10, it is also where it becomes less efficient than 
global reionization in suppressing the formation of satel- 
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Figure 8. Radial distribution of satellites in the fiducial 
fbk:sat/rei:G+L model. The observations include both the 11 
classi cal sate l lites and t hose detected in SPSS DR5 (compiled 
from iMated jl998f> and IKoposov et al.l J2009I) ), for which the 
SDSS DR5 threshold selection and sky coverage has been assumed 
(see text for details). The error bars illustrate the uncertainty due 
to Poisson statistics. Note that in this cumulative plot the errors 
are not independent. 

lites. Feedback alone (shown by dark blue dot-dashed lines) 
cannot explain the luminosity function in either model: the 
abundance of galaxies overshoots the data by a factor of 
~ 2 in the fbk:B06 case and by a factor of ~ 10 in the 
fbk:sat case. In both cases, supernova feedback is the domi- 
nant source of suppression of massive, bright satellites (with 
My < —10). As found in earlier studies, a combination 
of feedback and photoheating (global only in fbk:B06 and 
global + local in fbk:sat) is required to account for the abun- 
dance of low-mass satellites. 



4 TESTS OF THE FIDUCIAL MODEL 

We have selected the parameters of the fiducial model 
(fbk:sat/rei:G+L) on the basis of the luminosity function 
and Z — My of Milky Way satellites. In this section, we test 
the model against measurements of the radial distribution of 
satellites and the relation between the central mass density 
and luminosity. It is important to note that the model has 
not been adjusted during this comparison. 

4.1 The radial distribution of satellites 

Fig. [8] shows the radial distribution of the satellites that 
survive to the present day in the fiducial model for the six 
Aquarius halos. As outlined in Section [3] for magnitudes 
fainter than the classical dwarf regime (My > — 11), we 
select only those satellites that could have been observed in 



the SDSS DR5 survey. The detection of brighter satellites 
is assumed to be complete and so these are all included. 
Satellite distanc es and their un certainties (typically 10-20%) 
are taken from IMated (1199a) in the c ase of the classical 
dwarfs and from IKoposov et al.l (2009) in the case of the 
ultra-faint dwarfs detected in the SDSS DR5. 

Overall, the new model displays a radial distribution 
of surviving satellites similar to the observations (which 
have large uncertainties). The predicted distributions ap- 
pear slightly more concentrated than the data for 4 of the 6 
halos, less concentrated than the data for one and very close 
to the data for the remaining halo. 

4.2 The A/300 - L relation 

Dwarf galaxies are inefficient retainers of baryons and, as 
a result, are strongly dark matter dominated. This makes 
them ideal probes of the dark matter. In particular, their 
very central regions could contain informat ion about its 
identity (e.g. iNavarroT Frenk fc White] Il996h . Rough esti- 
mates suggest a mass of about 10 7 Mp within the visible 
parts of the Milky Way satellite (|Matedll99"a ; lGilmore et all 
l2007f l . More r obust analyses g i ve total masses of about 
3 x 10 9 M l|Wolf et all |2010| ; IWalker et all I2010I ). and, 
surprisingly, a common mass contained within the central 
300 pc, M300 ~ 10 7 Mq, independently of l uminosity over 
four orders of magnitude l|Strigari et al.ll2~00ct ). 

The high numerical resolution of the Aquarius halos al- 
lows us to calculate M300 directly from the simulations and, 
in conjunction with our modelling of galaxy formation in 
subhalos, to determine how it varies with satellite luminos- 
ity. Fig.[9]shows the M300 — L relation for the satellite galax- 
ies in the six simulated Aquarius halos, in the case of the 
fiducial model. We compute M300 directly from the simula- 
tions, by summing the masses of all particles within 300 pc 
of the centre of each subhalo. 

Since 300 p c is close to the resolut ion limit of the Aquar- 
ius simulations l|Springel et al.ll2008bl ), we have checked the 
numerical convergence of M300 using halo Aq-A, for which 
a higher resolution simulation, at level 1, is available. We 
created catalogues of subhalos matched in the level 1 and 
level 2 simulations by pairing up halos according to posi- 
tion, velocity and mass. We find an extremely good corre- 
lation between M300 for the most massive subhalos at the 
two levels of resolution, with more scatter present at lower 
masses, as expected. However, we find that the level 2 sub- 
halos systematically underestimate M300 relative to level 1 
by about 20%, even in the most massive systems. This dif- 
ference is sufficiently small that it will not affect our conclu- 
sions, but should be kept in mind when comparing model 
predictions with data. In the following, we present only the 
level 2 simulations. 

Fig. [9] shows that there is broad consistency in 
M300 between the A quarius simulations and the data of 
IStrigari et al.l (|2008l ). However, the masses in the simula- 
tions have larger scatter and, in some cases, slightly larger 
values than the observations. In all cases the simulations 
show a weak trend with luminosity that is not apparent in 
the data. The gasdyna mic simulations of ga l axy f ormation 
in Aquarius halo D by lOkamoto fc Frenk I (|2009h show a 
similar trend between M600 and luminosity. 

In some of our models, a population of very faint satel- 
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Figure 9. The M300 — L relation for satellites in the six Aquarius halos. Results for the fiducial model are shown as empty squares. 
Observational estimates of M300 (green stars) are from lStrigari et al,l i200Sl ). 



lites forms which have values of M300 lower than the typical 
~ 10 7 Mq estimated for the real satellites. It may be that 
systems like this exist in nature but have not been detected 
so far because of their intrinsically faint luminosities and 
low surface brightnesses. At the bright end the simulations 
show a large scatter in M300 above and below 1O 7 M0. Bet- 
ter observational statistics may help understand the reason 
for this discrepancy. 

The physical origin of the M300 - L rel ation has re- 
cently been discussed bv lStringer et al l (|2010h . We explore 



this question here by examining the dependence of M300 on 
satellite virial mass. 

4.3 The A/300 - Mvir relation 

In Fig. [TTJ] we plot the M300 — M v i r relation for satellites 
in the six Aquarius halos. Here, M v i r is the virial mass be- 
fore infall into the main halo. The solid line in each panel 
shows the M300 — M v i T relati on derived for dark matter halos 
with NFW density profiles l|Navarro. Frenk fc White Ill996l . 
Il997h . assuming the best-fit average mass - concentration 
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Figure 10. The A/300 — Mvii relation for satellites in the six Aquarius halos. Resu lts for the fiducial model are sho wn as empty squares. 
The solid line in each panel shows the M300 — M V - 1T rel ation determ i ned fr om the iNavarro, Frenk fc Whitel jl997T) dark matter profile, 
using the mean mass — concentration relation found bv lNeto et al~l 11200711 . 



relation found bv lNeto et ail (|2007l ). The M300 — M v i r re- 
lation has a similar shape to the M300 — L relation, that is, 
A/300 depends only weakly on halo virial mass. A power-law 
fit to the trend yields A/300 oc AfV 3 . This result can be un- 
derstood as follows. For a halo of mass A/ Vir = lO 9 A/ 111 
the ACDM cosmology, the ratio of 300 pc to r v \r is ~ 0.015 
and scales as 

/ \ -1/3 

15°° „ 0.015 [^r) . (4) 
r vir ^ lO 9 A/ J v ; 



At small radii (i n the limit of r <g; r B , where r B is 
the NFW scale radius; INavarro. Frenk fc White Ill997l ') in a 
given halo, the NFW density profile asymptotes to p oc r _1 , 
implying a mass profile M(< r) oc r 2 . Since the radius r3oo 
is much smaller than r s for the halo masses under considera- 
tion, the ratio A/300 /A/ V i r scales as (r3oo / V v ir) 2 . Substituting 
in the scaling for r3oo/Vvir in equation (4) yields the scaling 
A/300 oc A/V 3 found in the Aquarius simulations. 

The weak trend in the A/300 -A/ V i r relation is the under- 
lying cause of the weak trend in the M300-L relation and 



16 A. S. Font et al. 




Fi gure 11. Lumi nosity functions for the six Aquarius halos in 
the lLi et al] model. 



Figur e 12. Luminosity function for the Aq-A halo in the lLi et al] 
(2010) model with global reionization and feedback switched on 
and off. 



appears to be a robust prediction of cold dark matter the- 
ory unless baryonic processes in the forming dwarf galaxy 
are able to modify the inner regions of the dark matter den- 
sity profile significantly. 

4.4 A parameterization of the A/300 — L relation 

Summing up the results above, we can now understand why 
M300 appears to be independent of L, even if in reality it has 
some weak dependence. According to our results, L varies 



as M v 



where n ~ 2.5, and M300 varies as A/, 



1/3 



This 



yields the very weak dependence, M300 ~ lP , where /3 ~ 0.1, 
consistent with current observations. This trend is a robust 
prediction of our models which may, in principle, be tested 
by further photometric and spectroscopic observations of the 
faint satellites of the Local Group. 

Other theoretical studies find it similarly difficult to re- 
prod uce the completely flat M300 — L relation seen in the 



data dKoposoy et al|2009l;lMacci6 et al.|l2009l;lLi et al.ll2009l; 



Munoz et all 120091; lokamoto fc Frenk I 12003; iBusha et al 



2010l; Stringer et alj 20101 ') . In particular, iMunoz et al 



(2009) find a similar fit to ours (M 300 ~ L u 
-A/300 values derived from their simulation. 



") using the 



4.5 Comparison with the Li et al. (2010) model 

It is instructive to test the robustness of our model by 
comparing with the results of an independent semi-analytic 
model implemented in the same six Aquarius halos. This 
is the "Munich model" described by |Pe Lucia fc Blaizot] 
|2007t) and applied to ga l actic satellites by Li et al.h201fj ). 
The fee Lucia fc BlaizotJ (|2007l ) model reproduces a variety 
of observation al data in the local Unive rse and at high red- 
shift and, as |Pe Lucia" fc Helmil I 2008) have shown, with 



some modifications in the treatment of disc instabilities and 
star formation, it reproduces many observed physical prop- 
erties of our own Galaxy, including the age and metallic- 
ity distribution of stars in its different c omponents. (The 
modifications of Lucia fc Helmil (120081) do not alter the 
agreem ent with observations shown in |Pe Lucia fc BlaizotJ 
(120071) ). Additional updates of th e reionization a nd feedback 
prescriptions were introduced bv lLi et al.l l|2010l) to provide 
a better match to the observed properties of Milky Way 
satellites. A brief description of the relevant aspects of this 
model is as follows: 

- The cooling mo del used by | Li et al.l d2010l) is that 
originally proposed by IWhite fc Frenk ( 199ll ) and used in 
subsequent versions of the 'Munich model'. Specifically, the 
cooling rate depends on the temperature and metallicity 
of the hot gas, and is regulated by the infall rate in the 
'rapid accretion regime'. A comparison between this cooling 
m odel and that adopted by lBower et al.l (|2006l ) is discussed 
in |Pe Lucia et al.l (|2010l ) . Cooling via molecular hydrogen 
is not included, under the assumption that H2 is efficiently 
photodissociated. 

- The reionizatio n is modelled usin g the m e thodo logy 
of lCroton et~all (120061) . who ad opted the lbnedinl (|2000l ) for- 
malism. In the iLi et alj l|2010l ) model, however, the reion- 
ization is assumed to start at redshift zo = 15 and end at 
z r = 11 . 5, tha t is, earlier than implied by the calculation of 
Gnedin (2000), who obtaine d zp = 8 a nd z r = 7. The reion- 
ization scheme assumed in ILi et al.l (|2010l ) does not take 
into account the local photoionization. However, with the 
final choice of parameters, this model becomes similar to 
our preferred GALFORM model, fbk:sat/rei:G+L, which in- 
cludes both global and local reionization and which reionizes 
completely at z ~ 10. 
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Supernova feedback is modelled as in 
|Pe Lucia, Kauffmann fe White! (|2004T ). With the adopted 
parameters, galaxies with virial velocity Vvir < 87 km s _1 
have more h eated gas mas s per unit s tellar mass than in 
the st andard ICroton et alj j2006) and |Pe Lucia fc Blaizotj 
l|2007l ) models. Their feedback becomes inefficient for low 
mass galaxies, similar to the situation in our fbk:sat model 
(compare the blue dot-long dashed lines in Fig. [12] and in 
the right panel of F ig. 0). 

iLi et alj (|201dh applied this 'ejection' model to a series 
of high (but much lower than in Aquarius) resolution simu- 
lations of Milky Way-type halos. They found a good match 
to a variety of observations, including the satellite luminos- 
ity function, the luminosity-metallicity relation, the radial 
and size dis tribution s , and the central mass - luminosity re- 
lation. The ILi et al.l (|201Cf ) model yields somewhat higher 
metallicities than GALFORM for systems with Lv > 10 6 Lq 
(but is able to reproduce the metallicities of the ultra-faint 
satellites). 

We apply here the ILi et al.1 |2010l ) model to the six 
Aquarius halos analyzed in this study. Fig. [TT] shows the 
satellite luminosity functions, constructed in the same way 
as for GALFORM (see details in Section[3]). Within the scat- 
ter there is a reaso nably good agreemen t with the observa- 
tional estimates of iKoposov et all (|2008T ). 

Fig- HH shows the role of supernova feedback and reion- 
ization separately and is similar to Fig. [7] in the case of the 
GALFORM models. There are several common characteris- 
tics between this model and our fiducial model GALFORM 
model, fbk:sat/rei:G+L: firstly, feedback from supernovae 
is more effective at suppressing the formation of the more 
luminous systems, while reionization acts more effectively 
at the faint end (My ^ — 10); secondly, the reionization 
alone does not suppress galaxy formation sufficiently for the 
model to match the observed lu minosity function. (T his is 
reminiscent of the conclusions of lBenson et alj (|2002bl 'l who 
also found th at global reionization implemented using the 
iGnedinl (|2000h formalism has a relatively mild effect). 

The comparison between different semi-analytical mod- 
els highlights the existence of degeneracies in the way in 
which different physical processes - supernova feedback and 
reionization, in this case - affect the properties of the result- 
ing galaxies. However, as we have shown, it is often possible 
to break these degeneracies by comparing the model pre- 
dictions to a variety of observables rather than to a single 
property such as the satellite luminosity function. 



5 DISCUSSION AND CONCLUSIONS 

We have implemented a detailed treatment of reionization 
in the Durham semi-analytic model of galaxy formation, 
GALFORM, and applied it to the 6 high-resolution simu- 
lations of galactic halos of the Aquarius project. The UV 
flux produced by the galaxy population is calculated in 
a self-consistent way and the contribution from quasars 
is taken from the obs ervationally inferred spectrum of 
lHaardt fc Madaul l|l996l) . The UV flux inhibits star forma- 
tion by: (i) preventing gas accretion onto low mass halos, 
and (ii) offsetting the cooling rate of the gas already inside 
halos. These effects influence how galaxies form and evolve 
and this, in turn, affects the strength of the future UV back- 



ground, resulting in a self-consistent calculation of the cou- 
pled properties of the galaxy population and t he IG M. 

We use the formalism of lOkamoto et al~l ([2008 ) to cal- 
culate the suppression of gas accretion due to photoheating, 
and the CLOUDY software to calculate the net cooling rate 
of the gas (i.e., cooling plus heating rates) inside halos. We 
allow for an additional local UV flux (i.e., over and above 
that of the metagalactic UV background) generated by the 
progenitors of the Milky Way. This results in an earlier effec- 
tive redshift of reionization for the Milky Way and leads to 
further suppression of star formation in satellites, particu- 
larly in very low mass systems. However, in order to reionizc 
the Universe sufficiently early, our model requires an escape 
fraction of UV photons of nearly 100%. This is higher than 
current observational estimates. Although these are uncer- 
tain, it seems likely that the high efficiency of photoion- 
ization required by our models may reflect either the lim- 
itations of the inevitable approximations we have made to 
calculate the reionization process or the neglect of other pro- 
cesses such as supern o va-dr iven cosmic ray pressure which 
IWadepuhl fc Springell |201l[ ) argue could also play a role in 
suppressing galaxy formation in small halos. 

We find t hat a model with s upernova feedback as pa- 
rameterized bv lBower et al.l l|2006l ) predicts dwarf metallici- 
ties that are too low compared to observations. This can be 
rectified by assuming a lower efficiency of supernova feed- 
back in small systems. We have found that a model with a 
saturation in the supernova feedback efficiency (/3 = const) 
in dwarf satellite systems and with both local and global 
reionization provides the best match to the local luminosity 
function, the metallicity — luminosity relation, and the ra- 
dial distribution and central densities of present day dwarf 
satellites. As it retains the default feedback efficiency for 
more massive systems, this model also matches the large- 
scale luminosity function. 

In the fiducial model the suppression of faint satellites 
(My > —10) is achieved by a combination of supernova feed- 
back and (global+local) photoheating. However, the relative 
contribution of these processes varies across the range of 
dwarf galaxy sizes: supernova feedback is the dominant pro- 
cess for suppressing the more massive systems, photoheating 
plays the main role for the ultra-faint dwarfs. 

While the role of reionization in suppressing galaxy for- 
mation in small halos has long been recognized, the impor- 
tance of inh omogeneous reionization has only rece ntly been 
highlighted (|Busha et al.ll20ld ; iMufioz et al.ll2009f ). The lo- 
cal photoheating from Milky Way progenitors is a crucial 
ingredient in our model - unlike in most previous studies of 
the formation of satellites in the CDM cosmology, we find 
that global reionization by itself does not provide an ac- 
cep table solution to th e 'missing satellite problem'. While 
the I Bower et al.l (|2006l ) model gives a good match to the 
local luminosity function with global reionization only (see 
Fig. [4]), it does not match the Z — My relation and produces 
too few dwarf satellites when local sources of reionization 
are taken into account. In contrast, the saturated feedback 
model matches all of these datasets well. Global reionization 
occurs at z ~ 7.8 and locally at z ~ 10. 

In Appendix A, we show that reionization - with and 
without local photoheati ng - can be approxim ated using the 
standard v cut — z cut rule (Be nson et alj[2002bl ) . but with pa- 
rameters that differ from those typically adopted in other 
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semi-analytical codes. For example, in the fbk:sat/rei:G 
model, the best fit is given by the set (w C ut = 34 km s , 
Zcut = 7.8), and in the fbk:sat/rei:G+L model by the set 
(fcut = 34 km s _1 , z cut = 10). The z cu t values inferred 
this way are in very good agreement with the values derived 
from the full reionization treatment (see Fig. [3} and much 
less computationally intensive. 

A more detailed comparison between our semi- 
analytical model and other models c an be found i n Sec- 
tion|43]and Appendix B. The model of lLi et all (|2010h is the 
closest in methodology to the GALFORM model presented 
here and produces similar results. 

The recent renewal of interest in the formation of satel- 
lite galaxies seems to be leading to a consensus that it is 
possible to reproduce basic properties of the galaxy popu- 



tion and its effects on galaxy formation 


Busha et al] 


201C; 


Kodosov et al.l!2009; Maccio et al.l 201(1 


Busha et al.l 


2010; 


Li et al.ll2010f). Bovlan-Kolchin et al.l (201lt). however, have 



recently presented dynamical evidence that the largest sub- 
halos in ACDM N-body simulations are much too concen- 
trated to be able to host the brightest satellites of the Milky 
Way. It is unclear at present whether this discrepancy could 
be due to statistics l|Parrv et al.ll201ll ; dynamical data exist 
only for a handful of satellites in just one galaxy, the Milky 
Way) or whether it points to a more fundamental problem 
such as the existence of feedback processes not included in 
current models or even t o a different nature for the dark 
matter l|Lovell et al.ll201lh . 

In this work, we have asked the general question of 
whether, in the context of the ACDM cosmology, it is pos- 
sible to account for the observed properties of satellites not 
just plausibly but specifically within a broad model that 
has been developed to understand the galaxy population 
as a whole, not just the satellites. This is a more challeng- 
ing question than merely constructing a model restricted to 
satellites, but one that provides a stronger test of our un- 
derstanding of galaxy formation. We find, in agreement with 
earlier works, that a sufficiently high redshift of reionization, 
in combination with feedback from supernovae, can indeed 
reduce the number of satellites to a level which agrees with 
the data. The specificity of our model allows us to explore 
additional properties of satellites, such as their metallicity. 
We find that the metallicity- luminosity relation of satellites 
can only be explained at the same time as their luminos- 
ity function if supernova feedback saturates at low halo 
masses. Interestingly, our result raises the possibility that 
both the supernova feedback and reionization (global/local) 
are scale-dependent processes, which may have important 
consequences for our understandi ng of galaxy format i on. 

Like oth er recent models |Okamoto fc Frenk I |2009| ; 
iLi et al.ll2010T ). ours also provides an explanation for the ap- 
parent common mass s cale of dwarf galaxies . In agreement 
with the conclusions of IStringer et al.l (l2010h . we find that 
the weak dependence of M300 on L is a result of a weak de- 
pendence of M300 on M V i r . Our model predicts that there 
should, in fact, be a weak trend of increasing M300 with L 
and significant scatter in the relation at low luminosities. 
These predictions should be testable once improved mea- 
surements of M300 are available for larger samples of galax- 
ies. 

There have been significant and rapid advances over the 



past five years in both observations of dwarf satellites and 
our theoretical understanding of how they form. Despite this 
rapid progress, the conclusions of the previous generation of 
models remain essentially correct. The overall results seem 
independent of the details of any specific implementation of 
galaxy formation physics and lead to the conclusion that the 
broad visible properties of the satellite population are repro- 
ducible within the current cold dark matter paradigm. This 
represents an important success. Furthermore, the cold dark 
matter cosmogony makes clear predictions for the dark mat- 
ter content of dwarf satellites that are mostly independent 
of the bar yonic physics and are eminently testable by ob - 
servations jStrigari et al.ll2010l ; iBovlan-Kolchin et al.ll20llh . 
An important challenge that remains for galaxy formation 
theory is to verify if the good agreement with dwarf satel- 
lite properties can be retained while simultaneously match- 
ing the broader properties of much more massi ve galaxies 
and gal axies at higher redshifts (see, for example, iGuo et all 
l|2011al ) ) . Given the complexity of galaxy formation it is only 
by confronting models with such a broad range of data that 
a convincing theory of galaxy formation will emerge. 
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APPENDIX A: THE (Vcut,^cut) 
APPROXIMATION 

A self-consistent calculation of r eioniza tion in GALFORM 
was described by iBenson et al] (|2002al ). who showed that 
this calculation is well approximated by the (v cut , z cut ) 
model used subsequently in the GALFORM code. In this ap- 
proximation, the baryon fraction accreted by halos is given 
by: 

, _ f if z < z cut and v vir < « cut , , . 

Jb ~ \ (/ b ) otherwise, 1 ' 

where (/b) is the mean cosmic baryon fraction. This 'short- 
cut' greatly speeds up the semi-analytic code. It has been 
ado pted in several other semi-analy t ic mo dels su ch as those 
by iBullock. Kravtsov fc Weinberg! (|2000h and ISomervilld 
(2002). It is therefore interesting to check whether the new 
treatment of reionization implemented in this paper is still 
well approximated by a step function of the form above, 
and if so, to determine the best fit values of the (v cut , z cut ) 
parameters. 

We find that the simple (v cu t , z cut ) model still provides 
a good approximation to the full calculation of reionization 
even when local sources as included. We illustrate this with 
the fbk:sat model, although this conclusion is equally valid 
for other feedback schemes, including the default power-law 
feedback scheme, fbk:B06. 

Fig. [13] shows the (fcut, z C ut) approximation for 
the models with saturated feedback and self-consistent 
global reionization, with or without local reionization, 
fbk:sat/rei:G and fbk:sat/rei:G+L, assuming an escape frac- 
tion of ionizing photons of 100%. The luminosity functions 
are shown with different colours, each one corresponding to 
one of the six Aquarius halos. 

We find that the parameter values (v cu t, z cu t) = (34 km 
s _1 , 7.8) provide a good match to the model with global 
emission only, fbk:sat/rei:G. The best fit v cut is slightly 
higher than the critical circular velocity of 25 km s _1 ob- 
tained for the Okamoto, Gao, & Theuns (2008) accretion 
model. This is because in GALFORM photoheating not only 
prevents accretion onto halos below some mass but it also 
offsets radiative cooling losses (and therefore star forma- 
tion rates) for systems that are massive enough to accrete 
baryons). Note that z cut gives a very good approximation to 
the redshift of reionization calculated self-consistently (see 
Fig. 0). 

Adding the ionizing flux from the local sources turns 
out to be equivalent to shifting the redshift of reionization 
to earlier values, while keeping v cu t roughly constant. The 
best (Wcutj z cut ) approximation for the fbk:sat/rei:G+L case 
is (tw.Zcut) = (34 km s -1 , 10). 

We find that similar good (w C ut,Zcut) approximations 
exist for the other cases when / csc is varied. Perhaps more 
surprisingly, we find that the (v C ut,z cu t) approximation 
works well for a variety of other feedback schemes (includ- 
ing the default fbk:B06 and power-laws with varying slopes), 
regardless or whether or not a local photoionizing flux is in- 
cluded, and irrespective of the variety of merger histories 
represented in the Aquarius halos. It is also surprising that 
the approximation works quite well down to the regime of 
the ultra-faint dwarfs, even though it was first derived in 
the regime of classical dwarfs. The success of this approx- 



Milky Way satellites in ACDM cosmology 21 



1000 



_fbk:sat/rei:G (f =1) 

/ V esc ' 

,fbk:sat/(v =34,z =7.6) 



100 




imation could not have been predicted a priori, given the 
diffe rences betwee n our cu rrent models and those of the orig- 
inal iBenson et al.l (|2002aT ) model where the approximation 
was introduced. 



APPENDIX B: 
MODELS 



COMPARISON WITH OTHER 



Several other studies of Milky Way satellite galaxies in the 
context of the cold dark matter cosmogony have recently 
been carried out. The broad consensus from these inves- 
tigations is that theoretical models can achieve agreement 
with the observed distribution of satellite galaxy luminosi- 
ties given reasonable assumptions about the process of star 
formation and the ability of the reionization of the Universe 
to suppress galaxy formation in low mass halos. These stud- 
ies have included various pieces of physics thought to be 
relevant to the formation of low mass galaxies, but none 
have performed as detailed and physically complete calcu- 
lations as described in this work. Below we briefly outline 
the methods and results of these recent works and contrast 
them with our own. 
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Figure 13. Comparison between the models with a detailed 
calculation of reionization and the («outi*cut) approximation. 
Full lines represent fbk:sat/rei:G models in the top panel and 
fbk:sat/rei:G+L models in the bottom panel, respectively. The 
various colours correspond to the six Aquarius halos. Dotted 
lines in both top and bottom panels represent the models with 
the (i>cut , -Zcut ) approximation. The values for this pair of pa- 
rameters are: (v cu t , 2 cu t)=(34 km s _1 , 7.8) in the top panel and 
(«cuti ^cut)=(34 km s _1 , 10) in the bottom panel. 



5.1 Busha et al. (2010) 

iBusha et al.l (|2010l ) explored the effects of inhomogeneous 
reionization on the population of Milky Way satellites us- 
ing a combination of a high-resolution N-body simulation of 
the Milky Way halo and a lower resolution, larger volume 
simulation of the universe to assess spatial variations in the 
epoch of reionization. They find that reionization typically 
occurs between z — 6 and z — 12 for Milky Way-sized halos. 
To model the fo rmation of ga l axies within their dark matter 
only simulation, IBusha et alj (|2010| ) assume that halos must 
have reached a critical mass (corresponding to the atomic 
cooling limit of approximately 10 4 K) prior to reionization. 
They then assign luminosities to galaxies using one of two 
methods. In the first, the assume a one-to-one mapping be- 
tween luminosity and peak halo velocity, Umax- In the second, 
more physically motivated approach, a star formation rate 
proportional to a power of the halo mass is assigned to each 
halo between the time it reaches the critical mass and the 
epoch of reionization. The luminosity is found by applying a 
stellar population synthesis model and integrating over this 
st ar formation ra t e. 

IBusha et al] l|20ld) find that, with a suitable choice of 
star formation rate normalization, the second model pro- 
duces a good match to the faint end of the observed satellite 
luminosity function. (The bright end is underpredicted, but 
this is understandable as they do not allow for star forma- 
tion in massive halos after reionization.) Additionally, they 
find that their model produces a good match to the radial 
distribution of satellites while simultaneously producing val- 
ues of M300 in reasonable agreement with the data (although 
somewhat too high for the more luminous halos). 

In compariso n to our treatment, the model of 
IBusha et al.l (|20 lOh is much more simplistic. However, their 
results agree with ours at a qualitative level - sufficiently 
early reionization can reduce the number of satellites to a 
level compatible with observations. Additionally, their con- 
clusion that the abundance of satellites is very sensitive to 
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the epoch of reionization i s cons istent with our findings. For 
this reason, iBusha et al] (|2010h emphasize the importance 
of considering inhomogeneous reionization, consistent with 
our finding that localized ionization from progenitors of the 
Milky Way plays a key ro le in setting the local reionization 
epoch. IBusha et al.l (|2010h use the Via Lactea II simulation 
which has a halo mass of 2.6 x 1O 12 M0, somewhat more mas- 
sive than the halos considered here, but do not address the 
issue of whether their model produces a Milky Way galaxy 
with the correct properties (such as stellar mass). 

5.2 Munoz et al. (2009) 

iMunoz et aD (|2009h adopt a slightly more involved ap- 
proach, also utilizing the Via Lactea II N-body simulation. 
They consider four channels of star formation. In low mass 
halos at early times they allow stars to form via molecular 
hydrogen cooling. This process is stopped at z = 20 when 
molecular hydrogen is assumed to be dissociated. Prior to 
reionization at z ~ 11, stars are allowed to form in halos 
above the atomic cooling limit, while after reionization, star 
formation is restricted to higher mass halos. At late times, 
some further star formation resulting from metal line cool- 
ing is allowed. In each case, a simple parameterization is 
employed in which the mass of stars formed is proportional 
to the mass of the halo during the relevant epoch. Stellar 
population synthesis models are then used to infer galaxy 
luminosities. 

W ith a suitable choice of parameters, I Munoz et al.l 
(2009) find excellent agreement with the observed satellite 
luminosity function and show that molecular hydrogen cool- 
ing is important for producing the correct abundance of 
low luminosity satellites. They additionally explore the be- 
haviour of M300 as a function of luminosity and find broadly 
good agreement with the data although with significantly 
more scatter than observed and a trend for M300 to be too 
high in bright ga laxies and too low i n faint galaxies. 

The work of lMufioz et alj (|2009l ) contains many of the 
features of our own work (e.g., H2 cooling and the effect of 
photons generated locally). Our model, however, includes a 
rigorous treatment of the variety of physical process involved 
in galaxy formation such as gas cooling, star formation, feed- 
back and global reionization which have been tested against 
observations of the galaxy population as a whole at different 
redshifts. It i s uncl ear whether the ad hoc prescriptions of 
I Munoz et al.l ([2009 ) would lead to realistic galaxies beyond 
the satellites of the Milky Way. 

5.3 Koposov et al. (2009) 

The study bv lKoposov et al l (|2009h differs from other stud- 
ies in that it employs a semi-analytic method to follow the 
growth and evolution of the subhalo population inside a halo 
of final mass 1O 12 M . As we discussed above, it is possible 
that this mass may be too small for the Milky Way h alo. 
However, like most other studies, IKoposov et al.l (|2009h do 
not explore whether their model produces a central galaxy 
with a stellar mass comparable to that of the Milky Way. 
They explore a variety of prescriptions for assigning stars to 
subhalos, including simple models in which a constant frac- 
tion of baryons turns into stars and models in which stars 



can only form in halos above some critical characteristic ve- 
locity after the epoch of reionization at z — 11 (either with 
a sharp transition from star forming to non-star forming at 
this critical velocity , or a smoother transition motivated by 
the work of lGnedinl l|2000h ). 

They find that the faint end of the luminosity func- 
tion is made up of galaxies which formed their stars before 
the epoch of reionization, with brighter satellites forming in 
more massive halos after reionization. In fact, in their mod- 
els with a sharp transition at the critical velocity they see a 
bimodal luminosity function made up of these two types of 
galaxy. In models with a smooth transition the bimodality is 
lost and good agreement with the observed luminosity func- 
tion is obtained. Their results for M300 are consistent with 
those of other works (although with less scatter - a conse- 
quence of their neglect of scatter in the mass-concentration 
relation, as they note). 

The approach of IKoposov etail l|2009h is significantly 
more phenomenological than that described in this work, but 
it helps to confirm, as do other works, that an early epoch of 
reionization can plausibly suppress dwarf galaxy formation 
sufficiently to achieve agreement with the observations. 



5.4 Maccio et al. (2010) 

IMaccio et all (|2010h compare results from three different 
semi-analytic models of galaxy formation applied to high- 
resolution N-body simulations. Their halos span masses from 
1.2 to 3.6 x lO 12 A/ with a median of 1.7 x 1O 12 M , and so 
they could potentially suffer from the same shortcoming as 
the Aquarius halos of being of somewhat too low mass. Their 
semi-analytic models have been previously matched to the 
properties of the broader population of galaxies, but unfor- 
tunately, they do not specify whether or not they produce 
the correct mass of the central galaxy in these halos. 

In contrast to our approach, the subhalo informa- 
tion is not used to determine the evolution of satellite 
galaxies (e.g. to determine merging timescales). To add 
a suitable r eionization-indu ced suppression of galaxy for- 
mation, the iGnedinl {2000) filtering mass prescription is 
added to each model, with a reionization history taken from 
iKravtsov et alj (|2004r ) . They find that all three models can 
achieve a reasonable match to the observed satellite luminos- 
ity function with a reionization epoch of z — 7.5. However, 
they note that the original filtering mass prescription over- 
estimates the suppressing effects of reionization. Adopting 
the currently favoured suppression (which becomes effective 
in halos with characteristic velocities below ~ 30 km s _1 ), 
they find that a higher redshift, z = 11, of r eionization is 
requir ed to restore a good match to the data. IMaccio et al.l 
(2010) explore the roles of various physical ingredients in 
their models in achieving this match. In particular, and in 
agreement with this work, they find that the inclusion of 
supernova feedback is crucially important - without it far 
too many lumino us galaxies ar e form ed. 

The work of IMaccio etafl (|2010h is closest to our own 
in terms of the range of physics modelled and the detail of 
the treatment. However, these models lack the potentially 
important effects of molecular hydrogen cooling and do not 
include of a self-consistently computed reionization history. 
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5.5 Guo et al. (2010) 

IGuo et all (|2011aT ) use a semi-analytic model of galaxy for- 
mation to study the properties of galaxies in dark mat- 
ter halos spanning a wide range of scales by utilizing 
dark matter halo merger trees from the Millennium and 
Millennium-II N-body simulations. Subhalo information is 
taken from the simulations and used to track the merg- 
ing of satellite galaxies. Of relevance to this work, their 
Millennium-II merger trees resolve halos down to a mass of 
approximately 2 x 10 s Mq, which is sufficient to just resolve 
Milky Way dwarf satellites. To study Milky Way analogues, 
IGuo et all (|2011aT ) select halos from the Millennium-II sim- 
ulation which contain a disk-dominated (judged in terms 
of stellar mass) central galaxy with a total stellar mass in 
the range 4 x 1O 1O M < M 4 < 6 x 10 10 M o which results 
in a median halo mass of 1.44 x 1O 12 M . This once again 
raises the potential issue of the halo masses being too small, 
although in this case if further observational evidence sug- 
gests a higher mass for the Milky Way halo it would neces- 
sitate a recalibration of the IGuo et all l)2011al ) galaxy for- 
mation model to reduce the stellar mass of galaxies in dark 
matter halos of given mass. The galaxy formation models 
used are based on models previously used to successfully 
model vari o us aspe cts of the galaxy population. Importantly, 
IGuo et all (|2011aT ) explore how their model performs not 
only for the Milky Way satellite population, but also for the 
broader population of galaxies. They demonstrate that their 
model provides a good match to field and cluster galaxy lu- 
minosity functions while simultaneously matching that of 
Milky Way satellites. 

The key physics of reio nization is incorp orated into the 
model of IGuo et all ^Ollal ) in by use of the lGnedinl (|2000h 
filtering mass prescription, with a filtering m ass as a function 
of red s hift extracted from the simulations of lOkamoto et al~l 
( 20081 IGuo et all (|2011al ) find that they obtain good agree- 
ment with the luminosity function of Milky Way satellites 
using their standard reionization prescription correctly pre- 
dicts the number of bright satellites, but is marginally in- 
consistent with the number of faint satellites (producing 
somewhat too many, although the model remains plausi- 
ble given current uncertainties in the observational sample). 
They find that removing the effects of reionizati on only af- 
fects th e abundance of fainter galaxies. Although IGuo et al.l 
l|2011al ) do not discuss this point, this implies that super- 
novae feedback plays a major role in inhibiting the formation 
of satellite galaxies, in agreement with our own findings. 



